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Abstract -We have produced and observed metastable solid helium-4 below its melting pressure 
between 1.1 K and 1.4 K. This is achieved by an intense pressure wave carefully focused inside 
a crystal of known orientation. An accurate density map of the focal zone is provided by an 
optical interferometric technique. Depending on the sample, minimum density achieved at focus 
corresponds to pressures between 2 and 4 bar below the static melting pressure. Beyond, the 
crystal undergoes an unexpected instability much earlier than the predicted spinodal limit. This 
opens a novel opportunity to study this quantum crystal in an expanded metastable state and its 
stability limits. 
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■ Crystalline solids can be brought into metastable state 
^ with respect to fusion if surface melting is avoided. Ex- 
O perimentally this is achieved by heating small metal or 
I dielectric crystals embedded in stronger matrices beyond 
the bulk melting temperature Metals are the most 

frequently studied materials. Fast heating of bulk solids 



by subpicosecond laser pulses has also provided evidence 



J — . 6i such metastable states . Numerous computer sim- 
00 Illations have been used to unveil mechanisms limiting the 
^-H itietastability domain |21[S]-[5] . 

\^ • Solid helium offers novel opportunities. At low temper- 
^^atures (0 < T < 1.4 K), hep solid helium-4 only exists 
^^for pressures higher than P,„(T) ^ 25 bar. The melt- 
. . ipg line is nearly horizontal in the T-P plane, in contrast 
^ with other simple solids where it is nearly vertical. In this 
|v>< unique situation, pressure is the efficient control parame- 
5_] ter to explore the metastability region below the melting 
5^4ne. In addition, solid helium offers the peculiarity that 
quantum zero point motion plays a more important role 
than thermal motion. For melting, it does not obey the 
empirical Lindemann criterion [5]H Metastable expanded 
solid helium-4 has also been considered in the past in view 
of achieving the Andreev-Lifshitz-Thouless scenario for su- 
persolidity [111112] by lowering the vacancy energy [13] . Al- 
though this scenario is now considered less relevant [14115] , 



the interstitial and vacancy energies in a stretched lattice 
has been computed [14] , and a limit of the metastable re- 
gion has been proposed to lie 35 bar below the melting 
line [12]. None of these prospects ever came to an experi- 
mental realizatiorl^ and a presumably large domain for the 
solid metastable phase remains unexplored. 

In this letter, we provide full experimental evidence 
that metastable solid helium-4 can be obtained well be- 
low Pm(T) by a focused pressure wave inside bulk solid. 
This technique was used to study the liquid phase of he- 
lium in a metastable state with respect to its vapor or to 
the solid phase [TBlfT^ . For the hep lattice of solid hclium- 
4, the sound velocity is anisotropic. A pressure wave front 
originating from a point source is not a sphere, but a pro- 
late surface along the c— axis [201. Hence hemispherical 
transducers do not produce focused sound waves in that 
case. Decompression of hep solid helium with an hemi- 
spherical transducer has actually been tried with only a 
modest success [21] . Here we use an appropriate shaping of 
the ultrasound emitter (piezo-electric transducer PZT) to 
compensate for the sound velocity anisotropjU. The den- 
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0.26 at melting instead of 



^Attempts to obtain metastable solid helium at much higher pres- 
sures and temperatures using temperature as the control parameter 
are reported in 1171 

■^A cylinder of bulk ceramic PZT401 from Morgan Electroceram- 
ics was grinded on a CNC lathe prior to metallization and pol- 
ing. The inner surface matches half the wave surface of the quasi- 
longitudinal sound mode. The outer surface is parallel to the inner 
one providing a constant thickness of 2 mm. The equatorial part is 
then cut away over 1 mm in order to make the focal region visible. 
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sity modulation, below and above the coexistence density, 
is accurately measured using an interferometric imaging 
technique. 

CCD 




pulsed laser "k = 532 nm 



Fig. 1: A convergent pressure wave is produced in hep solid 
helium at 1.1 K by the excitation of a well shaped piezo- 
transducer (PZT). The density modulations near focus are 
monitored using a Jamin interferometer. The lens L images 
the focal region onto the CCD camera 

The experimental cell containing solid helium is cooled 
by a helium-4 cryostat with four optical ports. The tem- 
perature can be regulated from 1.1 K to 2.1 K. Internal 
pressure is measured by a capacitive gauge |22] calibrated 
along the liquid-solid coexistence line against an external 
pressure gaugcQ. The symmetry axis of the PZT is ver- 
tical and is taken as z-axis (Fig. [1]). The orientation of 
solid helium should be such that its symmetry axis (per- 
pendicular to the c-facet) coincide with the z-axis. He- 
lium density modulations are monitored via the refractive 
index using a Jamin interferometer (Fig. [I}. One of its 
arms passes through the acoustic focal region and is taken 
as the y-axis, while the other one crosses the cell in an un- 
perturbed zone. The focal region is imaged onto a CCD 
camera using an appropriate lens. The imaged surface 
is about 4 x 2.4 mm^. The measurement technique has 
been described elsewhere |23j and may be summarized as 
follows: optical phase shifts (^(2:, z) induced by density 
modulations are extracted from several images taken with 
different phases of the reference beam. This process is then 
reproduced for different values of the delay t between the 
sound pulse emission and the imaging laser pulse to con- 
struct a phase map movie. Because of the rotational in- 
variance of the sound field around the z-axis, the variation 
of the refractive index 5n depends only on the coordinate 
z and on the distance r to the z-axis. The optical phase 
map is then the integral along the y-axis of 5n{r, z), i.e. its 
Abel transform. An inverse Abel transform |24j of </'(a;, z) 
provides the radial profile 5n{r,z). Density maps can be 
computed from 5p(r,z) = po 5n{r, z) / {uq — 1) where po 
and no are the unperturbed density and refractive index 
of the solid. This method of determining sound fields has 
been quantitatively tested in water by comparison with 
an optical fiber hydrophone |23] . For an easier apprecia- 
tion of the modulation achieved, conversion to pressure is 

"Keller PAA-35/30bar/80797 



made using: 

6P[r,z)=x'-^ (1) 

where x = 269 bar [5S] is the bulk modulufH and no = 
1.036 [37]. This formula is the linear approximation of the 
equation of state, valid up to 5P c:i 5 bar with an error 
less than a few percent. 

When producing solid under depression, nucleation of 
the liquid phase on the PZT's surface must be avoided. Let 
5PpzT{t) be the acoustic pressure at the PZT surface, Tq 
the temperature of the cell and PmiTo) the static melting 
pressure. We note Pq the uniform pressure in the cell and 
(5Po ~ Po ~ Pm{To) its initial offset from the melting line 
(Fig. Po is chosen so that the condition SPpzrit) + 
6P0 > is always fulfilled, even at the troughs of the 
pressure wave. 
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Fig. 2: Experimental path in helium-4 phase diagram to ob- 
tain a crystal at To and Po- A ^ B: liquid pressurization at 
constant T . B: nucleation of an oriented crystal seed. B C: 
warming of the crystal seed in equilibrium with the liquid 
up to Tc ~ 1.4 K. C; growth of the crystal until the whole 
cell is filled. C — > D{To,Po)'- isochoric cooling. Coexistence 
curves are taken from Grilly [5S]. Solid dots are experimental 
data provided by the inner pressure gauge. In this example, 
5Po = 0.54 bar. 

To obtain an oriented and homogeneously pressurized 
crystal at (Tq, Po), a crystal seed with its c-facet horizon- 
tal is first grown at a temperature below the roughening 
transition of the c-facet (1.3 K) [25]. The horizontality of 
the facet is controlled optically in both x and y directions 
with 1° accuracy. Filling the cell with the crystal and pres- 
surizing it from Pm{To) to Pq at constant temperature Tq 
would be harmful for crystal quality. Indeed, because he- 
lium is a very compressible solid, material would have to 
be brought into the cell including into the PZT cup. This 
would result in plastic deformation of the crystal, possibly 
breaking it into several differently oriented pieces. Pres- 
sure gradients could also be created. Instead we follow the 
procedure shown in Fig. [5] The cell with the seed in equi- 
librium with the liquid is warmed up to some temperature 
Tc chosen so that P,n{Tc) > Po > Pni(Po)- This warm- 
ing along the melting line is done slowly enough to keep 
a quasi-static equilibrium between the pressure in the fill- 
ing line and the melting pressure at the cell temperature. 
The crystal is then grown at Tc until the cell is filled and 

^The use of formula JTJ is justified because solid helium com- 
pressibility is nearly isotropic (cn -|- C12 = C33 + C13). See 1261 
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the filling capillary is bloeked by the solid phase. Even- 
tually, the whole cell is slowly cooled at constant volume 
down to (Tq, Pq) in a little more than two hours. This slow 
isochoric cooling is expected to avoid crystal breaking, as 
well as temperature and pressure gradients. 
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Fig. 3: Refractive index map Sn{r, z) of a convergent pressure 
wave in licp helium single crystal. PZT voltage \s V = 185 V 
and propagation time from the PZT surface is t = 15.53 /is. 
For an easier appreciation of the wave field around the focus, 
the map is supplemented on the left of the axis by a mere 
symmetry. 

The experiment then runs as follows. A four-oscillation 
resonant RF pulse (frequency 1.25 MHz, amplitude V tun- 
able from 10 to about 250 V) is applied to the PZT and 
generates a convergent pressure wave train in the crys- 
tal. Phase maps are recorded for delays t ranging from 
10 /xs to 20 /.is with 0.05 /iS time steps. Radial refractive 
index maps 5n[r,z,t) are obtained from them. Such a 
map is shown in Fig. [31 From those maps, one can find 
the acoustic focus Zf as the location where 5n(0,z,i) has 
its maximum amplitude. In Fig. 21 the refractive index 
variation at acoustic focus (5n(0, is plotted versus 

time. The pressure scale of the right axis is computed 
using equation ([l}. The pressure reaches its minimum 
at train = 15.53 /.ts and its variation from the static 
pressure Pq is SPmin = —2.80 bar ±0.12 bar. The er- 
ror bar is computed from the r.m.s. noise on dn and 
corresponds to one standard deviation. The initial off- 
set pressure being SPq = 0.87 bar, one concludes that 
the solid is decompressed under its melting pressure by 
AP = SP,nin + SPa = -1.93 ± 0.12 bar. In another run 
with SPo = 0.23 bar, the value AP = -2.2 ± 0.05 bar was 
obtained. 

Does helium stay in the solid phase? If a liquid bubble 
had nucleated at the focus, it would have expanded for 
about a quarter of the sound period at a velocity close to 
the sound velocity. Thus it would have reached a diameter 
of the order of 100 /xm. Taking into account the refractive 
index difference between liquid and solid phases (about 
0.0035), the liquid bubble would have created a phase shift 
as large as 4 radians on a size larger than the 20 /lui im- 
age resolution (to be compared to maximum phase shifts 
of about 1 radian induced by the pressure wave). More- 




Fig. 4: Time evolution of refractive index and pressure varia- 
tions at focus for V — 185 V. Maximum transient depression 
{SPmin ~ —2.8 bar ±0.2 bar) is obtained at tmin = 15.53 us 
and is clearly below the static melting line indicated by the 
dotted hue at -SPo = -0.87 bar. 



over, a large asymmetry would have been found between 
positive and negative pressure swings. None of these two 
phenomena is observed. Thus transient melting docs not 
occur. 
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Fig. 5: Minimum depression achieved below the melting pres- 
sure AP versus the PZT driving voltage V. Diamonds and 
open circles correspond to different runs with different crystals 
and different initial pressure offsets SPo- Error bars correspond 
to one standard deviation of the noise on 5n. For the last point 
on each run, a defect appears at focus (see inset), which is also 
responsible for the larger noise. The crystal stability limit in 
AP is between the two dash-dotted lines. 

Is the crystal modified by the pressure wave? Recent ex- 
periments [29l[30] show that creation and motion of dislo- 
cations are involved in the response of solid helium-4 even 
to very low shear strain (e « 10~*). The temperature used 
here (1.1 K) is significantly higher than the temperature 
where the dislocations become free to move. Hence, it is 
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likely that with the strain produced in this experiment 
(typically 2x10-"* at the PZT's surface at 200 V, 4x10-3 
at focus) the solid becomes full of dislocations. These dis- 
locations can organize themselves into grain boundaries, 
leading to a subdivision of the crystal into slightly dis- 
oriented crystallites [21] ■ We looked for such a crystal 
modification in two ways. First, we have performed a set 
of experiments where the driving voltage {V) is gradu- 
ally increased from a reference voltage Kof = 60 V up to 
220 V. In between each different driving voltage acquisi- 
tion, a test acquisition at the reference voltage was taken. 
We then compared the successive test phase maps to the 
reference map and found them quite similar. For instance, 
the minimum depressions AP for the reference voltage are 
within 20% as shown in Fig. [5] (open circles). For compar- 
ison, a misaligned crystal with its c-axis at 20° from the 
PZT axis or a solid made of a few crystals give very dif- 
ferent density maps. As a second test, after sonication at 
100 V, the crystal was melted carefully at its coexistence 
pressure, reversing the procedure used for its growth. The 
liquid/solid interface reveals the existence of grain bound- 
aries around the focal region. The conclusion is that the 
local disorientation of the crystal is not sufficient to mod- 
ify significantly the pressure wave propagation. 

In the same series of measurements, we looked for the 
maximum achievable depression. It was found that the 
crystals break when V exceeds a threshold lying between 
185 V and 220 V. An extended defect with a different den- 
sity appears at the focal point (see inset of Fig. The 
corresponding value AP is between —1.93 and —2.40 bar. 
This threshold zone is determined as follows. For a given 
run at a given temperature, its upper limit is the lowest 
pressure achieved without breaking the crystal. To esti- 
mate its lower limit, we extrapolate AP from the last two 
points where the crystal remains stable up to the voltage 
where the crystal is broken. For the two runs of FiglSj 
the given threshold zone is the union of the corresponding 
threshold zones. 

Once defects are created in the crystal, they disappear 
in a few seconds after the sound pulse. However, subse- 
quent sound focusing is degraded, indicating a substan- 
tial crystal reorganization in the central region. Hence to 
carry out a new experiment, the crystal has to be melted 
again and a new one regrown using the procedure shown 
in Fig. [5J Clearly the limit encountered here is not the 
limit of the solid phase discussed in ref. [TB] which is ex- 
pected at much lower depression (—35 bar). In order to 
get some hints about its origin, 24 different runs were per- 
formed varying the temperature, the initial pressure offset 
(5Po, the repetition rate of sound pulses. Each run was 
made with a new crystal. It appears that AP somewhat 
depends on the crystal, although they are grown in similar 
conditions. Actually, AP may change by a factor 2 from 
one sample to another. It was found that AP is sensi- 
tive to the emission rate of sound pulses received by the 
crystal. A smaller rate tends to promote larger AP. Sur- 
prisingly, we did not find a strong correlation between 5Pq 



and AP. Finally, AP does not depend much on tempera- 
ture from 1.1 K to 1.4 K as shown in FiglHl It represents 
the lowest thresholds achieved for five temperatures in this 
temperature range. 
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Fig. 6: Threshold zones (crystal stability limit, see text) at 
different temperatures between 1.1 K and 1.4 K 

From this graph, it seems unlikely that the instability is 
related to the hcp/bcc coexistence line or its extrapolation 
in the metastable region. Note that these data were taken 
at slower sound pulse rates than the data of FiglS] This 
may explain the lower minimum pressures achieved. 

A further indication appears in FiglS) The values of AP 
do not scale linearly with V above F = 60 V. This reveals 
some non linear process in the production or the propaga- 
tion of the pressure wave, possibly dislocation creation or 
motion. The initial dislocation density may most probably 
vary from one crystal to another; their proliferation might 
depend on the number of strong sound pulses sent through 
the sample. If confirmed, this would appear as a technical 
limitation of the ultrasound technique to achieve larger 
depressions using a bowl shaped PZT. One can neverthe- 
less try to double AP by using two such PZTs, assembled 
face to face, to make a nearly complete wave surface. 

To summarize, we have demonstrated the production 
of metastable solid helium between 1.1 K and 1.4 K us- 
ing pressure as a control parameter. For amplitudes up 
to -2.0 bar below the melting line, the ultrasound tech- 
nique used does not seem to alter drastically the crystal. 
An unexpected instability of the solid phase was found. 
Nevertheless, a metastable region for hep solid helium-4 
is now accessible where crystal properties can be inves- 
tigated. Among these, bulk compressibility can be ad- 
dressed through Brillouin scattering. Vacancy energy will 
be more difficult to measure, but no less interesting. 
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